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Abstract 

Isospin violation effects in nuclei are investigated within a microscopic ap- 
proach which takes into account the inffuence of the residual strong interaction 
on the Coulomb interaction. The renormalization of the off-diagonal Coulomb 
field acting within a nucleus, by the residual strong interaction is calculated 
analytically in a simplified RPA. From the expression for the suppression coef- 
ficient of the isospin violating part of the Coulomb field, the isospin violating 
spreading widths of isobaric analog states are derived. The resulting reduction 
of the width is in agreement with the data on the isospin symmetry restoration 
and with some earlier evaluations of the isospin violating spreading widths. 
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The off-diagonal part of the Coulomb potential in nuclei plays an important role in deter- 
mining some properties of nuclear structure [|I|- |Tl|. Many effects of isospin mixing including 
the widths of isobaric analog resonances are governed by off-diagonal matrix elements of the 
Coulomb one-body potential. It was realized [Q]- [^] that the off-diagonal one-body Coulomb 
matrix elements are strongly renormalized by the strong interaction in nuclei. This renor- 
malization leads to reduced isospin mixing and the isospin symmetry is preserved to a better 
degree than what one would expect from unrenormalized Coulomb interactions. The method 
in which the renormalization was treated is by introducing collective effects into the radial 
motion in nuclei. The notion of the collective giant isovector monopole state (IVM) was 
introduced and its role as a mediator of many Coulomb mixing effects [Q- was discussed 
twenty five years ago. The reduction in the Coulomb matrix elements is achieved by intro- 
ducing a collective RPA upwards shift of the IVM and by a reduction of its strength due to 
the repulsive nature of the particle-hole interaction in the isovector channel. 

The problem of isospin symmetry violation in complex nuclear states attracted attention 
for many years [|T],^|TU|,|rT],^,^ in the general context of theory of chaotic systems The 



case of isospin symmetry plays a special role in this respect because the very microscopic 
origin of its violation is well known and therefore constitutes a good example where one can 
test the theories of chaos in many-body systems. 

A natural way to achieve this is to use an explicit realistic strong interaction that de- 
scribes well nuclear properties and to employ RPA to evaluate the one-body Coulomb matrix 
elements. We are able to develop approximate analytical formulas which describe the renor- 
malization of the off-diagonal Coulomb potential due to the strong nuclear force acting in 
nuclei. When realistic nuclear forces are used, the resulting Coulomb matrix elements are 
quenched and and are of the size that is capable to explain the experimentally observed 
reduction of the isospin violating spreading widths. 

We start with the nuclear Hamiltonian H written in the form 

H = H^ + Vs + Uc , (1) 
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where the first term Hq = p^/2m + Us{r) is the single-particle Hamiltonian of nucleons 
moving in a single-particle strong potential Us{r). We will use the harmonic oscillator form 
Us = ^^^^r^. The term Vs is the residual two-body strong interaction. The one-body Coulomb 
potential is denoted by Uc and we use the uniform density form: 

Uc = 3ZeV(2i?)(l -rV(3i?^)),r < R 

Uc = Ze^/r, r> R, (2) 

where R, a, and r are the nuclear radius, diffuseness, and the radial variable correspondingly. 

We will be interested below in the RPA renormalization of the off- diagonal part of the 
Coulomb potential, Uc, defined as the rest of Uc after subtraction of Uc which is diagonal 
in the representation in which the single-particle Hamiltonian Hq is diagonal 

Uc = Uc + Uc. (3) 



We now proceed |13;T3 by defining an anti-Hermitean operator A and perform the following 
unitary transformation. The transformed Hamiltonian is: 

H* = e^He-^ = 
Ho + Vs + Uc + Uc + [A, Ho] 

+ [IVs], (4) 

where we have used the decomposition and neglected all terms above the first order in 
the Coulomb interaction. To obtain the effective one-particle off-diagonal Coulomb field, U^, 
that results from the initial one, Uc, due to the resudual strong interaction effects, we should 
find the operator A in such a way that the single-particle off-diagonal Coulomb contribution 
in e^He~^ will be cancelled. The last term in Eq.(^) is a two-body operator. We introduce 
the decomposition: [A, V5] = ([A, V5])-|- : [A, V5] :, where the first is a single-particle term 
obtained by averaging over the nucleon occupied states, and the second two-body term, 
: [A, Vs] :, is the rest which yields zero under such averaging. 
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We can compensate the one-body Coulomb type field in the Hamiltonian H*, if we require 
the operator A to satisfy the equation 

Uc + [A, Ho] + {[A,Vs]) = 0. (5) 

Then, the transformed Hamiltonian takes the form 

H = Ho + Vs + Uc+:[A,Vs]: (6) 

where no single-particle off-diagonal Coulomb potential is present. In this case, the effective 
off-diagonal Coulomb potential field can be found from the equation 

U*c = -[A,H% (7) 

where the effects of the renormalization of the Uc are incorporated into the operator A [pT^ . 
In the language of the present approach, the operator A in Eq.(^ creates a small distortion 
of the nuclear density matrix when acting on the density matrix which contains no effects 
due to the residual interaction V^. Equations (|§),(0) can be seen to be equivalent to the main 
equation of the Theory of Finite Fermi Systems 

^ik = ^ik + ^iys)pa,ik- —^ik (8) 

a/3 £/3 - 5a - a; 

which relates an initial off-diagonal field to its effective field JF*, obtained via summation of 
the series of RPA diagrams involving the residual interaction and the particle- hole propagator. 
A |T^. In Eq.(P), indices refer to the single-particle nucleon states which label the matrix 



elements of the fields JF, JF*, the interaction Vs and the single-particle occupation probabilities 
Hi of the states with the energies e^. In our case, the external (Coulomb) field is static, so 
uj = Indeed, substituting Eq.(|^) into Eq.(|) and collecting the terms one obtains an 
equation analogous to Eq.(|^) with JF* and JF replaced by U* and U. Equations (|§),(0) 
provide a convenient way to obtain analytical results for the renormalization of the potential 
in the approximation which we develop below. 
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Up to now, we did not specify the explicit form of the residual strong interaction. To 
obtain results for the renormalization according to Eqs.(|^,0), one needs to specify V5. We 
use below the Landau-Migdal parametrization [|T^, |T6[ of the strong interaction V5, 

Vsir,v')=C6iv-r')[f + f'TT' + gaa' + g'TT'aal (9) 

which has been widely used for heavy nuclei (see ||15|, ||16|). One therefore considers the form 
(P) as a reliable parametrization of that reproduces correctly the main properties of actual 
residual interaction. The strengh constants f ,f',g, and g' depend on the nuclear density. The 
relevant part of the interaction Eq.(^, Vr, that is responsible for the renormalization of the 
one-body off-diagonal Coulomb field is given by the first two terms of Eq.(P) whose explicit 
form reads as 

Vrir, r') = C6{r - r') - (/e. - /.„)(p(r) - po)Po ' + 

HfL - ifL - fDipir) - Po)Po'Krr')]. (10) 

Here, p(r) is the nuclear density, and po denotes its value in the center of nucleus. Here 
C = 300 MeV fm? is the universal Migdal constant. The values of the strength parameters 
are [|Tg: /e. = -1.95, /i„ = -0.075 f',, = 0.05 fl^ = 0.675. 

From symmetry considerations, the operator A can be taken in the form: 

A = XpiiVpfp) + (fpVp)) + x„((V„f„) + (f„Vn)), (11) 

to match equation (^) with some constants Xp, Xn to be determined. 
In the expression for {[A, Vs]) we have, after some calculations. 




The off-diagonal part of Eq.(|T2[) which is to be substituted into Eq.(|) is given by the off- 
diagonal parts of the second and the fourth terms in ( [T^ ) involving the off-diagonal terms of 
Tp and (denoted by and f,^). The following combinations are introduced: 
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The new constants are defined as follows 
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Now, substituting Eq.(|T2D into Eq.(|) and accouning for (|l5),(|l5) we obtain, after separating 
the similar operator structures, the following system of equations: 

^ 2A\Uo\\2-^PP i^pp) 2A|l/o|V2 P" l^pn) 

_z_Cpo(5p -^F') \-IL£po(^F -'yF' 

2A|C/o|V2 "P l-^np) ^ 2yl |i7o| V2 ™ ^ nr 

The solution for this simple system of linear equations is found inverting the matrix in Eq.([Tf 
the value of the determinant is given by the formula 
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Instead of the initial off-diagonal Coulomb field 
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using Eqs.([7D,(PTD and finding Xp, x„ from the system of equation (|T^), we obtain now proton 
and neutron effective (renormalized) fields which are given by the expressions: 
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(18) 
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(the asterisk marks the renormahzed quantities). The effective field that acts on neutrons is 
non zero, it contains now the contributions of the proton-neutron component of the residual 
strong interaction. The sign of this contribution is opposite to that of the proton field. 

For the renormalization of the isovector part of the off-diagonal Coulomb field, [/*'■*•' = 
Ij^ip) _ IJ^'^\ that can be written as 

we obtain from Eqs.(0,0) the following expression: 

I _q*p-q*n _ 



S q 



2 _ ^ \2 

pn) 



(19) 



This analytical result gives the value of the suppression factor S of the isospin violating part of 
the Coulomb field as a function of the nuclear charge and mass, and the interaction constants. 
A simple numerical evaluation of Eq. (^) for the case ~ Z and the values of the interaction 
strengths presented above (see Eq.(0)) gives the following value of the suppression factor S 

5 =1.8 ±0.4 (20) 

where the uncertainty is mainly due to that in the quantity 7. The obtained value of the 
reduction factor squared 5^ ~ 3.4±1.4 can be compared to the suppression factor that results 
from comparison of the experimental spreading widths of the IAS resonances with the values 
calculated using the unrenormalized Coulomb potential: 

Sip ~ 3 - 4 (21) 

Therefore, the present theoretical result for the renormalization of the off-diagonal isovector 
Coulomb field can be considered as a microscopic theoretical explanation of reduced isospin 
symmetry breaking in nuclei. 

The utilization of the above effective off-diagonal Coulomb potential is limited to the 
inside of the nucleus and to the states localized inside. We should keep in mind that the 



reduction derived above is applicable only to the part of Eq.(0) and is not applicable for 
the outside (r > R) part that goes as K 

The way in which our present result could be applied is illustrated in the following two 
examples: 

A. Isospin mixing in ground states. In order to calculate isospin mixing we can use the 
potential U* in the expression [Q: 

. .^2 ^ |(n/j|?7*W|n + l/j)p 

where \nlj) are single-paticle states and \n + llj) the corresponding radially excited s.p. 
states. The denominator is the energy difference between the unperturbed single-particle 
states. The collective effects due to RPA correlations are taken into account in . We find 
that 

where e is the amount of isospin mixing in the single-particle model without the suppression 
effect calculated in the present work. We see that our model predicts a reduction in isospin 
mixing of about a factor S*^ ~ 4 compared to the pure single-particle model in which the 
renormalization of the off-diagonal Coulomb field due to the strong interaction is not taken 
into acount. 

Ji. Spreading width of isobaric analog resonances. In the single-particle approximation the 
spreading width is given by 

= -Im{A\U,Gt,Uc\A) 

where 1^4) is the analog state, and Gp is the optical potential model Green's function. As 
already discussed this kind of calculation overstimates the spreading widths of isobaric analog 
resonances by factors close to 5 for standard optical potentials [|I|. If we replace Uc by f/*^*-* 
we find that the spreading widths are much closer to the experimental ones and close to the 
results obtained in Ref. p where the collectivity of the IVM is taken into account. 
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We should stress that this reduction apphes only to the spreading widths and not to 
the decay (escape) widths r\. The latter involves continuum wave-function which are not 
localized inside the nucleus and therefore T\ depends on the 1/r part of Uc- As mentioned 
above the reduction we find does not apply to 1/r and therefore it affects the escape width 
much less, in agreement with the calculations in Ref. 

Some of the reduction of isospin mixing occurs already in the Hartee-Fock (HF) approxi- 
mation 0] . The attractive proton- neutron force couples the proton and neutron distributions 
reducing the difference in the two caused by the Coulomb force. In an unrestricted HF calcu- 
lation that allows for deformations and for charge-exchange modes one obtaines in the wave 
fuctions components that can be expressed as a particles coupled to the various RPA states 
(including the IVM) of the core. Thus in such extended HF calculation effects of the RPA will 
be included. However the use of such unrestricted RPA is not practical and the use of such 
scheme for the quantitites calculated in this work is limited for the following reasons, (a) In 
order to incorporate correctly isovector effects in the HF one must include charge-exchange 
excitations, which is a difficult task, (b) The HF calculation in an odd A system involves 
additional approximations, (c) The usual HF approach in nuclei with a neutron excess intro- 
duces spurious isospin mixing and the symmetry potential acts as an isospin breaking term 
1^. (d) The HF approximation cannot be used for unbound states. 

To summarize, we have developed here an approximate analytical approach to study 
the role played by the residual strong interaction in reducing the Coulomb matrix elements 
violating isospin symmetry in nuclear states. We calculated for the first time analytically 
the RPA renormalization of the off-diagonal Coulomb potential starting from the nuclear 
Hamiltonian with the residual two-body strong nuclear forces. When realistic nuclear forces 
(Landau-Migdal parametrization) is used our approach produces quenched Coulomb matrix 
elements that are capable of explaining the experimental isospin violating spreading widths. 

This work is supported by the Grant for Basic Research of the Israeli Academy of Science. 
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